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TECHNICAL NOTE 3616 


CHARTS FOR ESTIMATING ROTOR-BLADE FLAPPING MOTION 
OF HIGH-PERFORMANCE HELICOPTERS 
By Robert J. Tapscott and Alfred Gessow 

SUMMARY 


Theoretically derived charts shoizing steady-state and first- and 
second-harmonic flapping coefficients are presented for helicopter 
rotors having hinged rectangular blades with a linear twist of 0°, 

-8 , and -l6°. The charts^ showing the rotor flapping coefficients for 
combinations of inflow ratio and blade pitch angle, are presented for 
tip-speed ratios ranging from O.O5 to O.5O. 

The theory on which the charts are based does not contain any 
small-angle assumptions regarding blade-section inflow angles and veloc- 
ities. These charts are therefore considered applicable, mthin the 
limits discussed, for flight conditions involving high inflow velocities 
and large regions of reversed velocity. Other than including an approx- 
imate allowance for stall in the reversed-velocity region, the charts do 
not include stall or compressibility effects. 


INTRODUCTION 


A knowledge of the steady-state flapping behavior of lifting rotors 
is useful in the design of helicopter hubs and control systems, in esti- 
mating rotor-fuselage clearances, and as a prerequisite to the numerical 
(that is, step-by-step integration) evaluation of the aerodynamic charac- 
teristics of rotors. Equations for calculating rotor-blade flapping 
are available from various theories; however, in order to reduce compu- 
tations to a minimum, this paper presents charts whereby theoretical 
flapping values can be obtained directly. 

The charts were prepared from the extended rotor theory of refer- 
ence 1. Since the equations of reference 1 do not assume the rotor 
inflow angle or the rotor angle of attack to be small, they can be used 
in calculating rotor-blade flapping over a iTide range of operating con- 
ditions. The assumption is made, howevfer, that, outside the reversed- 
velocity region, the section angles of attack are small; therefore, these 
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angles can be replaced by their sines. Also, -vrLth the exception of an 
approximate allowance for stall in the reversed-velocity region, the 
theory and, thus , the charts of this paper do not include stall or com- 
pressibility effects. The charts were calculated for hinged blades 
having 0°, -8°, and -l6° of twist and for a range of conventional values 
of blade mass constant. 


SYMBOLS 


P blade flapping angle at particular azimuth position 

measured from plane perpendicular to axis of no 
feathering, radians or deg 

aQ constant term in Fourier series that expressed p> 

hence, rotor coning angle 

aj^ coefficient of cos nijf in Fourier series that 

expresses p 

bj^ coefficient of sin ni]f in Fourier series that 

expresses p 

7 mass constant of rotor blade j expresses ratio of air 

forces to mass forces cpaR^/lj^ 

p tip-speed ratio, ^ 

A inflow ratio, 7 


difference betireen blade-root and blade-tip pitch 
angles, positive when tip angle is larger, deg 

Qq blade-section pitch angle at blade root, deg 

blade-section pitch angle at O.75 radius; angle 
between line of zero lift of blade section and 
plane perpendicular to axis of no feathering, deg 

tt/ blade-element angle of attack at tip of retreating 

' ’ ‘ ^ blade at 270° azimuth position, deg 


“(u^=0.4)(270°) 


blade-element angle of attack at radius at which 
tangential velocity equals 0.4 tip speed at 270° 
azimuth position, deg 
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V 

a 


°r 

V 

r 


a 

R 

P 

a 


velocity of helicopter along flight path, fps 

rotor angle of attack; angle between axis of no 
feathering (that is, axis about which there is no 
cyclic-pitch change) and plane perpendicular to 
flight path; positive when axis is inclined rear- 
ward, deg 

blade-element angle of attack, measured from line of 
zero lift, radians 

induced velocity at rotor (always positive), fps 

radial distance from center of rotation to blade ele- 
ment, ft 

rotor angular velocity, radians/ sec 

blade radius measured from center of rotation, ft 

mass density of air, slugs/cu ft 

slope of ciirve of section lift coefficient against 
section angle of attack in radians (asstimed equal 
herein to 5* 73) 

mass moment of inertia of blade about flapping hinge, 
slug-ft^ 


If 


blade azimuth angle measured from downwind position 
in direction of rotation, deg 


c ' blade-section chord, ft 

Uj, nondimens ional component of resultant velocity at 

blade element 


ME?THOD OF MAYISIS 


The analysis of reference 1 was developed for linearly twisted 
rectangular blades mth the flapping hinge located on the rotor center 
and perpendicular to the rotor axis and to the blade span. All veloci- 
ties, angles, forces, and moments are referred to the axis of no feathering. 
An explanation of this system of axes, together with means for applying 
the equations based on this system to pure feathering or combination 
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flapping -feathering configurations, is given in reference 2. Configu- 
rations involving offset or tilting of the flapping hinges can he eval- 
uated hy the numerical method presented in reference 5^ which can also 
he used to compute flapping motion in flight conditions outside the range 
of the analytical theory. 

As used herein, hlade flapping is the angle hetvreen the blade-span 
axis and a plane perpendicular to the axis of no feathering. Since, in 
the theory, the hlade is considered inflexible, the calculated flapping 
angle should he considered to represent the angle between the plane of 
no feathering and a line through the flapping hinge and the 
station of the blade. The flapping angle is expressed as a function of 
the azimuth angle hy means of a Fourier series expansion that can be 
terminated for most pirrposes after the second harmonic. The expression 
for the flapping angle is then 

3 = a^ - a^ cos ijr - b^ sin ^ - a^ cos 2i]r - bg sin 2\]f 

The coefficients are related in reference 1 by means of five simulta- 
neous equations expressed as functions of tip-speed ratio p, inflow 
ratio A, pitch angle 0q, blade tvrlst and airfoil-section charac- 

teristics. The specific values of these variables that were used in 
the present paper in preparing the charts for rotor flapping coefficients 
correspond to the values represented in the performance charts of refer- 
ences 4 and 5» the five simultaneous equations were solved for the same 
range of variables. 


DESCRIPTION OF CHARTS 


The rotor coning angle aQ and the first-harmonic flapping coef- 
ficients a^ and bj^ are presented in figure 1 as a function of A 
and for tip-speed ratios ranging from O.05 to O.50 and for blade 

tvTists of 0°, -8°, and -l6°. In the figure, Uq and b^ are divided 

by the blade mass constant 7; thus, all the coefficients are essentially 
independent of 7. 

As in any theoretical calculation that does not make allowance for 
blade stall, it is desirable to include limit lines on the plots to 
indicate the flight conditions at which the theory is expected to become 
inaccinrate as a result of stall and, at the same time, to indicate 
limiting conditions of operation. Such limit lines are included in the ' 
perfoimance charts of references and 5 ^ foi" example, and can be used 
in conjunction with the flapping charts of this paper. For the sake of 
convenience, however, the inboard and outboard limit lines of references 
4 and 5 are included in figure 1. Also, to avoid crowding of the lines. 
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the limit lines for each of the twist values are shoira on separate plots; 
the limit lines apply to 0° twist on the Uq// plots, to -8° twist on the 

plots, and to -l6*^ twist on the hQ_/7 plots. The limit lines corre- 
sponding to a particular twist that appear on one plot can he applied to 
the other plots at a given tip-speed ratio by simply matching the corre- 
sponding values of A and 6. 75* 

The second-harmonic flapping coefficients are presented in figure 2 
for values of tip-speed ratio of 0.30, O.ilO, and 0.50 for 0]_ = -8°. 

Values of the coefficients for tip-speed ratios less than 0.30 are not 
given because the magnitude of the coefficients for the lower tip-speed 
ratios is less than half a degree, which is considered to be iriLthin the 
accuracy of the theory. Also, second-harmonic flapping coefficients for 
0° and -16° of twist are omitted because they differed by a negligible 
amount from the values shorn for -8°. Because the blade mass constant 
cannot be factored out for the second-harmonic flapping coefficients as 
it can for the first-harmonic coefficients, figure 2 contains separate 
plots for 7 = 15, 7=9, and 7=3. 

The charts of figures 1 and 2 can be used directly to yield flap- 
ping coefficients- if A, 0.Y5-> Mown at the start. If, 

as more often is the case, the known starting quantities are tip-speed 
ratio, thrust coefficient, and power, then the corresponding values of 
pitch angle and inflow ratio can be obtained from the charts of refer- 
ences if and 5. 

It will be noted that the effects of twist on blade flapping are 
small and that the differences shown are probably within the acctiracy 
of the calculations. 


APPLICABILITY OF CHARTS 


The usefiilness of the charts in different applications wherein an 
estimation of blade flapping is required is, of course, dependent on the 
assun^jtions and limitations of the theory on which the charts are based. 
These limitations and the extent to which they affect the usefulness of 
the charts are discussed in the following sections. 


Assumptions and Limitations 

The extended theory from. which the charts presented herein were 
derived (ref. l) does not contain any small-angle assun^jtions regarding 
blade-section inflow angles and velocities. The assumption is made, 
however, that outside the reversed-velocity region the section angle of 
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attack Oj, is small and can "be replaced "by sin o^. Even with = 20 °) 
this assumption causes only about a 2-percent difference in the numerical 
value of the angle. For these reasons, it is believed that the charts 
are applicable for cases wherein the rotor inflow velocity is relatively 
large, such as might occur with high-performance helicopters flying at 
high foiward speeds or at hi^ rates of climb. 

It has been demonstrated (see ref. 2 , for exan5>le) that, for flight 
conditions which do not involve appreciable stall or conpressibility 
effects, analytical rotor theory similar to that used in preparing fig- 
ures 1 and 2 can predict the coning angle and the longitudinal and lateral 
flapping angles within about 1 °. In order to achieve this accuracy in 
predicting b^, however, it was necessary to asstune a longitudinal var- 
iation, instead of a uniform distribution, of induced velocity over the 
rotor disk. Inasmuch as the theory on which the present charts are based 
assumes uniform inflow, it is expected that the accuracy of the values 
of b^ given in figure 1 wo^lld be improved by applying a correction for 
nonuniform inflow. The correction to b^^ for nonuniform inflow given 
in reference 6 woiild be sufficiently accurate for most purposes. 

The charts presented herein were calculated for rotors having uni- 
form chord blades and zero flapping-hinge offset. In addition, although 
the coefficients a^ and b^ are divided by the mass constant, a fixed 
value was used in the calculation of these coefficients. This procedure 
does not eliminate the higher order effects of mass constant on these 
coefficients (aQ and bj_); however, subsequent checks indicated these 

higher order effects to be well within the accuracy of reading the charts. 
Thus, these charts, as were previous chart presentations of the theory 
of reference 1, are considered applicable for the conventional range of 
blade mass constants (approximately 0 to 25). 


Satisfactory conparisons are indicated also between calculated blade 
motion and blade motion measured for blades having taper ratios (root 
chord divided by tip chord) as hi^ as 3^1» Iii such cases, however, the 
taper ratio was accounted for theoretically by basing the rotor solidity 
and blade mass constant on an equivalent blade chord. For the determi- 
nation of rotor solidity, the equivalent blade chord is defined as 


2 

/ cr dr 

J_0 . 

So 


'R 


; whereas, for the determination of mass constant. 


0 


cr^dr 


r^dr 


r 

Jo 


r 5 dr 


is considered to be more accurate. 


The effect of flapping-hinge offset on blade motion appears to vary 
with both mass constant and tip-speed ratio. Calculations (ref. 
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example) show that these effects can he considerable. In the event that 
calculations are desired for cases with significant amounts of flapping- 
hinge offset (greater than 1 or 2 percent), the theory of reference 5 
may he used. In this regard, the charts of this paper should prove useful 
in providing initial estimates for the iterative process, and thus con- 
siderably reduce the time required to obtain a solution. 


Applicability to Seesaw Rotors 

Although calculated for rotors having individually hinged blades, 
the charts presented herein are also useful in the case of seesaw rotors; 
for example, the longitudinal flapping coefficient a^ is unaffected. 

With regard to the coning angle a^, and thus the lateral coefficient 

bj^, the validity of the value obtained from the charts presented herein 

will be dependent on the stiffness of the blade-hub combination relative 
to the stiffness imparted by the centrifugal forces due to rotation. 

For typical seesaw-rotor systems, the blade-hub stiffness appears to be 
about 8 percent, or less, of the stiffness imparted by the centrifugal 
forces; thus, the charts should give fair approximations for most pur- 
poses. In any case, the charts may be used to estimate the preset 
coning angle aQ for the seesaw-rotor system. It should be noted again 

in this connection that the charts as derived represent the angle bet^reen 
the no-feathering plane and a line through the flapping hinge and the 
station rather than a blade-root angle. 


Stall and Compressibility Effects 

As mentioned previously, the theory on which the charts of this 
paper are based does not make any allowance for blade stall or the var- 
iation of blade-section characteristics with Mach number. When the sec- 
tion angles of attack exceed the stall angle or when the Mach number for 
lift and drag divergence is exceeded, the charts are not expected to 
predict flapping coefficients accurately. Inasmuch as conditions for 
which representative blade-section angles of attack reach 12° and l6° 
are indicated on the charts, an estimate can be made of stall conditions 
of the rotor for given flight conditions. It is believed that the effects 
of moderate amounts of stall (tip angle at ij; = 270 ° exceeding the stall 
angle by as much as ll-°) can be estimated by the trends indicated in ref- 
erence 2. In cases where larger amounts of stall are present on the 
rotor or where a more extensive study of the effects of stall is desired, 
the method of reference 3 may he applied. 

Unpublished calculations indicate that the effect of compressibility 
on blade motion is to bring about a general increase in the values of the 
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flapping coefficients. The largest increment is seen to occur in the 
longitudinal flapping coefficient a^^; the increment increases ’tTith Mach 
number to about 1° at a Mach number of 1 at the tip of the advancing 
blade. The changes in coning angle a.Q and lateral flapping coefficient 

bj^ are about e4ual, and amount to about one-half of that noted for 

a^. Such trends, although generally small, may become significant for 
some purposes and should be considered when the charts of this paper are 
used to predict blade flapping for cases in which compressibility effects 
may be present. In such cases, the effects of compressibility may be 
evaluated quantitatively by use of suitable airfoil characteristics in 
conjunction with the step-by-step method discussed in reference 3* 


CONCIUDIWG REMAEKS 


Theoretically derived charts showing steady-state and first- and 
second-harmonic flapping coefficients have been presented for helicopter 
rotors having blades mth a linear tmst of 0°, -8°, and -l6°. The 
charts, showing rotor flapping coefficients for combinations of inflow 
ratio and blade pitch angle, were presented for tip-speed ratios ranging 
from 0.05 to 0.50. Inasmuch as the theory on which the charts are based 
does not contain any small-angle assumptions regarding blade-section 
inflow angles and velocities > the charts are considered applicable, 
iri-thin the limits discussed, for flight conditions involving high inflow 
velocities and large regions of reversed velocity. Other than including 
an approximate aULowance for stall in the reversed-velocity region, the 
charts” do not include stall or coDipressibility effects. 


Langley Aeronautical Laboratory, 

National Advisoiy Committee for Aeronautics, 
Langley Field, Va., October 51^ 1955* 
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Figure 1.- Rotor coning i 
as a function of pitcl 
-8°, and -l6° of twis' 
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Inflow ratio, X 


Flapping lines 
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Inflow ratio, X 
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Limit lines 

“(l.0)(270“)‘*^° 
“(ux=0.4)(270')^'^° 

.16 “(l.0)(270°)"'^° 

°(uy=0.4)(270')‘*®° 



0 I I I I I L \ \ ^ 1/ I I I I 

-.40 -.32 -.24 -.16 -.08 0 .08 .16 

Inflow ratio, X 

(b) n = 0.10. 


Figure 1.- Continued 
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Flapping lines 

0|=O° 

0,=-8° 

0|=-I6° 


Limit lines 

“(1.0) (270°)"'^° 

“(u-]-=0.4)(270“)'''^° 

“(1.0) (270'’)”’®° 

“(Uy=0.4)(270'’)"’®° 


Figure 1 .- Continued. 
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-.40 -32 -.24 -.16 -.08 0 .08 .16 

Inflow ratio, X 

Flapping lines 

0, = O° 

0,=-8° 

0,= -16° 

Limit lines 

“(1.0)(270°)^'^° 
“(u-r=0.4)(270°)'‘’^ 

“( 1 . 0 ) ( 270 °)='®° 

“(uj=0.4)(270°)'’^ 


Inflow ratio, X 

1.2 
1.0 
.8 

^,deg .6 
7 .4 

.2 
0 

-.40 -.32 -24 -.16 -.08 0 .08 .16 

Inflow ratio, X 

(f) \i = o,4o. 





Figure 1.- Continued. 


MCA TW 3616 



-.16 -.08 O .08 .16 

Inflow ratio, X 


i 

i 

■ 

B 

■ 

c 

r 


□ 

[I 


z 

z 


1 

■ 

m 

■ 

2 

0 “ 

\lote 
! 

1 

: Lim 

it Itn 

es Q 

pply 

to - 


ist. 

■ 

■ 


■ 

m 


■ 

1 







K 

■ 

■ 

m 

m 

■ 

i 

i 

B 








■ 

m 

m 

m 

B 

B 

a 

B 








B 

■ 

m 

B 

1^ 


m 

B 




■ 

■ 

■ 

■ 

m 

B 

■ 

B 

. 

m 

B 

■ 



■ 

■ 

■ 

■ 


E 

m 

B 

m 

B 

m 

B 

M 


■ 

■ 




Hi 

B 

Ei 

E 

Hi 

m 

B 

i 




■ 

Hi 



/ 





/ 





■1 

■ 





Z 



7^ 





■ 

■ 







/ 






■1 

■ 






7 






Rapping lines 


Limit lines 

~ “(tO) (270”) ' '2° 

“(uj"0.4)(270°)°I2‘ 

°(l.0)(270')''®° 

°(uy»0.4)(270'’)''®' 
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Inflow ratio, X 
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Inflow ratio, X 


(g) (1 = 0 . 50 . 


Figure 1 .- Concluded. 
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:36 -.32 -.28 -.24 -.20 -.16 -.12 -.08 -.04 0 .04 .08 

.y:;l 5 Inflow ratio, X 


ba.deg 0 


-.36 -.32 -28 -.24 -.20 -.16 -.12 -08 -.04 0 .04 .08 .12 

y=9 Inflow ratio, X 


I I I I I I I I I I I I I I I 

All values were within limits Oand -0.1 and therefore are omitted — 
I I I I I I I I I I I I I I I I r I I I 


-.36 -.32 -.28 -24 -.20 _.I6 -.12 -.08 -.04 0 .04 .08 .12 .l( 

y=3 Inflow ratio, X 

.(a) ji. = 0.30. 

Figure 2 .- Second-harmonic flapping coefficients as a fimction of pitch 
angle and inflow ratio for hlades having - 8 ° twist. 
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Qz, deg 
bg, deg 

Qa , deg 
ba, deg 

Qa, deg 
ba, deg 
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I i I I I I I I I I I I I I I I I I I I I 
All values were within limits 0 and-0.1 and therefore are omitted. 
J 1 I I 1 I I I I I I 1 I \ — I — L 


-.36 -.32 -.28 -.24 -.20 -.16 -.12 -.08 -.04 0 .04 .08 .12 

^3 Inflow ratio, X 


(•b) II = o.ho. 


Figure 2.- Continued. 











